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Abstract - Modeling results are presented for proton-induced 
degradation of charge-coupled devices (CCDs) used in the 
Advanced CCD Imaging Spectrometer instrument on the 
Chandra X-Ray Observatory. A methodology is described that 
provides insights regarding degradation mechanisms and on- 
orbit performance for front-illuminated and back-illuminated 
CCDs. Proton-induced changes in charge transfer inefficiency 
are modeled. The observed amount of on-orbit degradation can 
be accounted for using a proton spectrum at the CCD location 
that is reduced in magnitude by a factor of ~1E5 compared to the 
spectrum incident on the spacecraft. 

I. INTRODUCTION AND BACKGROUND 

T HIS paper presents results of modeling activities that 
examined proton-induced degradation of charge-coupled 
devices used in the Advanced CCD Imaging Spectrometer 
(ACIS) instrument on the Chandra X-Ray Observatory, which 
was launched on July 23, 1999. ACIS consists of an imaging 
array that includes four front-illuminated (FI) CCDs and a 
grating (or spectrometer) array that includes six CCDs, four 
of which are front-illuminated and two of which are back- 
illuminated (BI). Front-illuminated CCDs were observed to 
exhibit charge transfer inefficiency (CTI) degradation during 
the first month of ACIS operation. In contrast, back- 
illuminated ACIS CCDs did not degrade. Several workers 
subsequently described the Chandra radiation environment, 
the observed CCD degradation, modeling activities, and 
mitigation approaches [1]-[7J. The general conclusion 
reached was that a fraction of the protons entering the x-ray 
telescope aperture are scattered to the CCD location and 
cause performance degradation. Related studies were 
conducted for the XMM-Newton X-ray Observatory [7] and 
the ASCA x-ray astronomy satellite [8]. 

A detailed modeling methodology is presented herein that 
offers insights regarding degradation mechanisms and on- 
orbit performance for BI and FI devices used in the ACIS 
instrument. Proton-induced changes in charge transfer 
inefficiency are considered. The CTI modeling results are 
employed to provide an estimate for that fraction of the 
proton fluence incident on the spacecraft that reaches the 
CCDs. Our modeling of ACIS CCD degradation was 
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performed a few months after the Chandra launch [9J. The 
present work describes that modeling effort and then 
compares the results with those obtained by other workers. 

II. Device Structure 

Figures 1 and 2 present cross-sections for FI and BI CCDs, 
respectively, based on information compiled from several 
sources [10]-[13]. The OBF (Optical Blocking Filter) used in 
ACIS is also included in each cross-section since it provides 
some shielding of the CCDs to incident protons. The OBF 
dimensions shown are those for the filter used with the grating 
array of six CCDs. (For the four imaging-array CCDs, the 
thickness of the bottom A1 layer in the OBF is 0.13 pm. 
Thicknesses of the other two layers in the OBF were the same 
for the grating and imaging arrays.) Figure 3 presents the 
doping profile for the n-layer in FI and BI devices [1 1]. This 
profile yields a junction depth of -0.48 pm and a substrate 
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doping density of 2.7E12 cm , or 5000 Q-cm. Figure 4 
shows the simulated potential in the n-layer under the 
assumption of no fixed oxide charge and for a gate voltage of 
0 V [11]. This simulation provides an estimate of the depth 
into the n-layer of the middle of the buried channel (-0.32 
pm). A buried channel (BC) width of -0.1 pm is assumed. 

m. Modeling Methodology 

The proton energy spectrum incident on ACIS CCDs was 
uncertain when our modeling activities [9] were conducted. 
Work was performed subsequently to estimate that 
information [3], [7], In general, the spectrum incident on the 
spacecraft will be altered during transit to the CCDs through 
the effects of shielding and scattering. The proton fluence will 
be reduced substantially at the CCD location compared to the 
fluence incident on the spacecraft. The simplifying 
assumption was made in our modeling that the shape of the 
proton energy spectrum incident on the OBF is the same as 
that incident on the spacecraft. (To simplify the analysis, 
normal incidence on the OBF and the CCDs is also assumed.) 
It is further assumed initially that the fluence at a given 
energy incident on the OBF is also the same as that striking 
the spacecraft. We then deduce both the fraction of the 
external proton fluence and the effective acceptance angle for 
incident protons that would produce the CCD degradation 
observed in space. 
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Fig. 1. Cross-section of FI CCD. The OBF dimensions are for the filter used 
with the ACIS grating array. 
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Fig. 2. Cross-section of BI CCD. 


The modeling methodology involves transporting the 
proton spectrum through the filter and the intervening CCD 
structure until protons, with degraded energy, reach the buried 
channel in FI and BI devices. The average dose deposited in 
nonionizing processes (i.e„ the displacement damage dose, 

Drf) in the buried channel for both types of CCDs is then 
determined through detailed analyses. That dose is used in 
conjunction with space and ground measurements of CTI 
degradation to estimate that fraction of the proton fluence 
incident on the spacecraft that reaches the ACIS CCDs and 
the corresponding effective acceptance angle. Shielding 
effects of the HETG (High Energy Transmission Grating) 
used with ACIS are also considered. 



Fig. 3. Doping profile versus depth into the n-layer [1 1]. 
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Fig. 4. Potential profile versus depth into the n-layer [11]. 

IV. Proton Spectra 

Figure 5 shows the proton spectrum incident on the 
spacecraft plotted in terms of differential fluence in units of 
protons per cm 2 -steradian-orbit-MeV. This input spectrum 
was calculated using the Space Radiation code, which 
restricts the lower-energy limit to 100 keV. (The validity of 
the underlying AP8 model is questionable at lower energies.) 
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Transport calculations require knowledge of the range-energy 
relation, R(E), for protons incident on a silicon target. The 
SRIM code [14] was used to determine R(E), which is shown 
in Figure 6. To simplify the modeling, it was then assumed 
that all protons lose their energy in a manner consistent with 
the range-energy relation (i.e., no straggling). Thus, for a 
proton with energy E to reach a specific depth x in silicon, the 
criterion R(E) > x must be met. Upon reaching depth x, the 
proton will have an energy E x such that R(E) = R(EJ + x. In 
this manner, the input proton energy spectrum can be 
transported to any depth of interest. 
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Figure 5. Proton differential fluence spectra. The spectrum incident on the 
spacecraft is shown as well as that spectrum transported to the middle of the 
buried channel for FI and BI CCDs. The simplifying initial assumption was 
made that the proton energy spectrum and fluences incident on the OBF are 
the same as that incident on the spacecraft. 


It is important to determine the transported proton energy 
spectrum for energies well below 100 keV in device active 
regions. The nonionizing energy loss rate (NIEL), which 
correlates with the amount of displacement damage produced, 
increases with decreasing proton energy [15] (Figure 6), and 
thus low-energy protons contribute significantly to device 
degradation. This subject is discussed in [16], [17]. (Silicon 
NIEL values were available in the present study for proton 
energies of 1 keV and higher. Our calculations of Si NIEL for 
lower energies (100 eV to 1 keV) indicate that displacement 
damage doses reported below for the buried channel may be 
underestimates by ~10%.) Figure 5 also shows the differential 
proton spectra for FI and BI devices at the center of the 
buried channel, i.e., after transporting the input spectrum to 
that location in both cases. 

For FI CCDs, in our modeling approach 100-keV protons 
do not reach the buried channel because their range in Si- 
equivalent material is -0.92 |im, whereas the edge of the BC 
is at a Si-equivalent depth of ~1.3 pm (including the OBF and 
the top layers of the CCD). However, most of the protons 


with energy >100 keV reach the BC and deposit nonionizing 
dose therein. In the case of the BI CCDs, the distance 
required for a proton to reach the buried channel is tens of 
micrometers. Therefore, protons in the input spectrum with 
energy below 100 keV do not matter in the analysis of CTI 
degradation. (Protons reaching the buried channel in both FI 
and BI devices with degraded energies <100 keV all have 



Fig. 6. Rate of nonionizing energy loss [15] and range versus incident proton 
energy for silicon. 

V. nonionizing Energy Deposition 

The following approach was used to determine nonionizing 
energy deposition in CCD active regions. The input proton 
spectrum (Figure 5) was transported through the OBF and the 
various device layers (Figures 1 and 2) in the manner noted 
above. At each depth increment into the device, a proton 
energy spectrum was obtained, such as those shown in Figure 
5 for the middle of the buried channels in FI and BI devices. 
At each of those depth increments, the displacement damage 
dose was calculated by integrating the product of the fluence 
at that depth and the NIEL over the entire energy spectrum. 
That process was repeated through the entire active region of 
the CCDs. Figure 7 shows displacement damage dose versus 
depth determined in that manner. The curve shown is the 
nonionizing energy deposition profile for silicon. To use this 
curve for the other materials in the OBF and the CCD, an 
equivalent thickness of silicon must be determined based on 
the density of a given material relative to Si. The results 
presented in Figure 7 were used to obtain the nonionizing 
energy deposition profile in and near the buried channel. 
Figures 8 and 9 show the results for FI and BI CCDs, 

respectively. D^ varies by -10% over the buried channel 
width in the front-illuminated device, whereas for the BI CCD 
the variation is < 1%. 
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VI. Charge Transfer Inefficiency Degradation 

Dale et al. [18] demonstrated that changes in charge 
transfer efficiency (CTE) for proton-irradiated CCDs scale 
with NIEL. It seems reasonable to expect that the radiation- 
induced increase in CTI (termed ACTI), where CTI = 1- CTE, 
in any buried-channel CCD will scale linearly with the 
displacement damage dose deposited in the device active 
region, i.e., in the buried channel. Nonionizing energy loss in 
that region produces traps that degrade the transfer efficiency. 
It is well established that there is a linear dependence of ACTI 
on irradiating particle fluence (<(>) for both monoenergetic 
particles and a spectrum of particle energies. (For example, 

see [19], [20] and references therein.) Since Da is linearly 
proportional to <|> for a specific particle type and energy, the 
linear dependence of ACTI on Da follows. (In general, to 
analyze the effects on CTI of a spectrum of irradiating 
particle energies, one must determine the contribution to Da 
in the device active region at each energy.) The constant of 
proportionality between ACTI and Da is expressed here by 
Ken* 

ACTI = KcnDa (1) 
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Fig. 7. Displacement damage dose deposited in FI and BI CCDs. 


For the CCDs used on ACIS, ground measurements of CTI 
degradation following 40-MeV proton irradiation were made 
prior to the Chandra launch. The value for Ken obtained in 
that work was expressed [21] as (11 ± 2)E~5 per nonionizing 
rad, where one nonionizing rad equals 6.25E7 MeV/g. 
Converting to the present units yields a value for Ken of (1.8 

± 0.3)E-12 (MeV/g) . 40-MeV protons are quite penetrating 
and lose, relatively little energy in traversing the CCD. Thus, 
the NIEL value is well defined and equals that at 40 MeV 

2 

(i.e., 4.3E-3 MeV-cm /g). Consequently, the displacement 
damage dose deposited in the device active region is known 


accurately, resulting in an accurate Ken determination at that 
energy. 



Fig.8. Displacement damage dose deposited in and near the buried channel 
in FI CCDs in the grating array on ACIS. 



Fig. 9. Displacement damage dose deposited in and near the buried channel 
in BI CCDs in the grating array on ACIS. 

Low-energy proton irradiations of ACIS CCDs at 104, 157, 
and 390 keV were conducted at Goddard Space Flight Center 
using a van de Graaf accelerator [21]. The Ken value 
obtained for 390-keV proton bombardment is (8.4 ± 0.6)E-5 
per nonionizing rad [21], which converts to (1.3 ± 0.1)E-12 

(MeV/g) . The present modeling approach was used in 
conjunction with additional information in [21] to reassess 
that value. This reassessment was performed because of the 
possibility of displacement damage dose uncertainty. An 
incident 390-keV proton loses a significant portion of its 
energy before it reaches the buried channel region (-20%). 
Accurate determination of Ken requires careful determination 
of the displacement damage dose deposited in the buried 
channel. The present work yields a NIEL value of -0.19 
MeV-cm 2 /g for degraded 390-keV protons that reach the 
buried channel. The measured value for ACTI per unit fluence 
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is 4.2E-13 cm 2 at 390 keV [21]. Using the presently 
determined displacement damage dose per unit fluence (i.e., 
NIEL) with that ACTI value yields a revised value for Ken of 

(2.2 ± 0.2)E-12 (MeV/g) 1 at that energy. This approach 
yields good agreement between the Ken values for 40-MeV 
and 390-keV proton bombardment of ACIS CCDs (i.e., 1.8E- 
12 and 2.2E-12, respectively). In the analysis that follows, we 

employ Ken = 2E-12 (MeV/g) 1 for all proton energies. In 
general, Ken depends on the CTI measurement conditions 
(e.g., background charge, clocking conditions, temperature). 
We assume here that this Ken value is applicable for the on- 
orbit measurement conditions for ACIS CCDs. 

Applying the present modeling approach to determine Ken 
at 157 keV based on CTI measurements given in [21] yields a 

value on the order of IE- 12 (MeV/g) V We consider this 
value to be in reasonably good agreement with that obtained 
at 390 keV and 40 MeV in view of the significant analytical 
uncertainties in the 157-keV case. In our analysis, a 157-keV 
proton loses nearly all its energy (-90%) before reaching the 
buried channel, and exhibits significant further degradation as 
it traverses that channel. Thus, NIEL values in the buried 
channel and the resulting displacement damage dose 
deposited therein are very dependent on the device 
dimensions assumed and on other modeling assumptions 
(normal incidence, no straggling, no spectral variation, etc.). 

Consequently, the error in the calculated Dj, and thus in Ken, 
is relatively large for incident 157-keV protons. 

In the present one-dimensional analysis, based strictly on 
the range-energy relationship 104-keV protons do not quite 
have enough energy to reach the buried channel in FI devices, 
even with the OBF removed, so no degradation is predicted. 
However, CTI degradation was observed at 104 keV in 
ground irradiations [21] performed with no OBF in place. 
This discrepancy between modeling predictions and 
experimental results may be due to the various assumptions 
made here. Those assumptions are critical for proton energies 
where the particle range is comparable to the location of the 
device active region. To explore this apparent discrepancy, 
we performed Monte Carlo transport calculations for 104-keV 
protons using SRIM [14]. Those calculations demonstrated 
that when straggling is accounted for, very-low-energy 
protons do reach the buried channel and will cause CTI 
degradation. However, we were unable to quantify the 
displacement damage dose deposited because accurate NIEL 
values are not known for those energies. Another possible 
reason for the apparent discrepancy is proton energy 
variations for the irradiation source used. In the present 
analysis, we assumed that the incident protons used in ground 
experiments were monoenergetic. To achieve a low-energy 
proton beam using a van de Graaf, a scattering foil is often 
used to degrade a higher energy beam. If the input beam 
energy is large compared to the desired output energy, the 
spread in the output energy spectrum can be significant. This 
possible spread in proton energies incident on the CCD would 


need to be accounted for carefully in a modeling approach 
such as that described herein. 

VII. RESULTS AND DISCUSSION 

When the ACIS CCDs were exposed to the space radiation 
environment for eight Chandra orbits, the pre-irradiation CTI 
of roughly IE-6 increased to -2E-4 in FI devices [1], [13]. 
Using Eq. (1) and the Kcri value determined above of 2E-12 
(MeV/g)' 1 , a displacement damage dose of ~1E8 MeV/g 
deposited in the buried channel would cause that amount of 

degradation. We now compare that value with the D<j 
calculated herein for the FI CCD buried channel, which 
assumed that the proton spectrum and fluence incident on the 
spacecraft is also incident on the OBF and CCDs. That 
comparison then allows an estimate to be made of that 
fraction of the incident proton fluence that reaches the CCDs. 
From Figure 8, the calculated displacement damage dose at 
the middle of the FI CCD buried channel is -9.4E10 MeV/g- 
sr-orbit. For eight orbits, a dose of -7.5E11 MeV/g-sr is 
obtained. For the dose deposited in the BC to be 1E8 MeV/g, 
as determined above, then the effective solid angle that admits 
protons to the CCD region needs to be -1.3E-4 sr. Under this 
interpretation, the fraction of the omnidirectional proton 
fluence that reaches the CCD is then ~(1.3E-4)/4rc, or-lE-5. 
The proton spectrum is assumed to remain unchanged, but the 
fluence is reduced by that factor. 

For BI devices, based on Figure 9 the displacement damage 
dose deposited in the buried channel during eight orbits is 
-1.1E9 MeV/g-sr, as compared to -7.5E11 MeV/g-sr for FI 
CCDs. Applying the solid acceptance angle calculated above 

(1.3E-4 sr), we obtain an estimated Dj for the BI buried 
channel of -1.5E5 MeV/g. From Eq. (1), that dose is 
predicted to increase the CTI by 3E-7, which is negligible. 
Many orbits (-5000) would be needed to deposit sufficient 
dose to cause the CTI in BI CCDs to increase to 2E-4. 

The effect of the HETG is now considered. After the FI 
CCDs were observed to experience a CTI increase to -2E-4, 
there were two additional orbits in which the HETG was 
inserted and the HRMA (High Resolution Mirror Assembly) 
covers were open. For this situation, the HETG provides 
additional shielding of the CCD beyond that already 
considered. Using the same methodology as employed above, 
we calculated the effects of that shielding on the displacement 
damage dose deposited in the buried channel of FI CCDs. 

The result is a D<j per orbit of ~9E6 MeV/g. This is 9% of the 
dose needed to cause the CTI to increase to 2E-4 (i.e., 9% of 
1E8 MeV/g). Thus, after an orbit with the HETG in place the 
CTI is predicted to increase by 1.8E-5. The resulting total 
CTI, -2.2E-4, most likely would be noted as a negligible 
change compared to the previous measurement of 2E-4, in 
agreement with results given in [13]. (Note that proton 
scattering in the HETG will decrease the fractional CTI 
contribution to < 9%.) 
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Kolodziejczak et al. [3] also analyzed proton-induced 
degradation of the ACIS FI CCDs. The differential proton 
spectrum they used was based on AP8MAX, as compared to 
AP8MIN in the present study. (In the critical energy range 
between 100 and 500 keV, the AP8MIN fluence values we 
used are -50% higher than those obtained using AP8MAX 
and employed in [3].) Kolodziejczak performed scattering 
analyses to determine the solid acceptance angle for incident 
protons to reach the ACIS CCDs. For the case of simply 
scaling the input proton spectrum (i.e., fluence-scaling only) 
to achieve the observed CTI degradation, they obtained an 
acceptance angle of 1.6E-4 sr, which agrees well with the 
value of 1.3E-4 sr obtained here. That close numerical 
agreement may be somewhat fortuitous in view of the 
differences between the analytical methodologies employed 
and the input proton spectra. To first order, this agreement 
shows that one can account for ACIS CCD degradation based 
on the damage produced by incident protons that reach the 
focal plane through a very narrow energy-independent 
effective acceptance angle. 

In the scattering analyses performed by Kolodziejczak et 
al. [3], they determined the acceptance angle as a function of 
proton energy. When that dependence is accounted for, their 
modeling underestimated the observed CTI degradation by a 
factor of 3 to 4. That apparent discrepancy led them to 
suggest that other processes may be important, such as 
damage produced by heavier ions. They also note that such 
conclusions depend on the accuracy of their scattering 
analyses. 

Nartallo et al. [7] also analyzed proton scattering for the 
Chandra x-ray telescope. For a fixed acceptance angle, over 
the proton energy range from 100 keV to 1.5 MeV, they 
obtained -9E-6 as that fraction of the incident proton fluence 
that reaches the CCDs. (The energy units on Figure 3b in [7] 
are evidently mislabeled as keV instead of MeV.) Their result 
agrees closely with the value of -IE-5 obtained herein. This 
agreement indicates that the present simplified modeling 
approach provides a physically reasonable description of 
ACIS CCD degradation. 

VIII. CONCLUDING REMARKS 

A methodology for determining proton-induced 
degradation of FI CCDs used in the ACIS instrument on 
Chandra is presented in this paper. Results of CTI analyses 
and modeling are described. The results are consistent with 
on-orbit observations of ACIS CCD degradation. In summary, 
charge transfer inefficiency degradation in CCDs used in the 
ACIS instrument on Chandra appears to be$ accounted for 
quite well quantitatively in terms of nonionizing energy 
deposition in the buried channel by relatively low energy 
protons. 
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